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Abstract: The frequent glitches (sudden increases of the apparent an- 
gular velocity) observed in certain pulsars are generally believed to be at- 
tributable to discontinuous angular momentum transfer to the outer neutron 
star crust from a differentially rotating superfluid layer, but the precise mech- 
anism is not quite elucidated. Most explanations invoke vortex pinning as 
the essential mechanism responsible for the build up of strain in the crust 
that is relaxed, either by fracture of the solid structure or by discontinous 
unpinning, during the glitch. It is shown here that there is another mech- 
anism that could give rise to strain, and subsequent fracture, of the solid 
crust, even if vortex pinning is ineffective: this is the effective force arising 
from the deficit of centrifugal buoyancy that will be present whenever there 
is differential rotation. This centrifugal buoyancy deficit force will be com- 
parable in order of magnitude, but opposite in direction, to the force that 
would arise from vortex pinning if it were effective. 
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1 Introduction. 



The ultimate motivation for this article is the problem of explaining one of 
the salient observational features of isolated (non-binary) pulsars, which is 
that comparatively long periods of continuous "spin down" of the observed 
frequency Q are occasionally interrupted by small "glitches". Such a glitch 
consists of a sudden small increase, Sfl say, that partially cancels the contin- 
uous negative variation Afl that has been accumulated since the preceding 
glitch. 

Since very soon after its discovery in 1968, it has been generally agreed 
that the pulsar phenomenon is attributable to a strong magnetic field an- 
chored in the outer crust layers of a central neutron star. The observed 
frequency Q is to be interpreted as the rotation frequency of the outer crust 
layer, whose continuous spin down is evidently due to the continuous de- 
crease of the angular momentum J due to radiation from the external mag- 
netosphere. After thirty years of work, two basic problems remain. 

The first is to account for the spectrum (from radio to X-ray and beyond) 
and the detailed pulse structure of the radiation, which are presumed to 
depend on the still very poorly understood workings of the magnetosphere. 

The second problem - the one with which the present article is concerned 
- is to account for the frequency "glitches". It is generally recognised that 
the glitches must be explained in terms of what goes on in the interior of the 
neutron star, and it is also generally believed that the glitch phenomenon is 
essentially related to the property of solidity that is predicted (on the basis 
of simple, generally accepted theoretical considerations) to characterise the 
crust of the neutron star after it has fallen below the relevant extremely high 
melting temperature, which occurs very soon after its formation. 

The purpose of this article is to draw attention to the potential impor- 
tance, as a mechanism for glitches, of the stresses induced in the crust just 
by the effective force arising from the deficit of centrifugal buoyancy that will 
be present whenever there is differential rotation. 

It is to be noticed that centrifugal buoyancy is a phenomenom that has 
been previously considered in the context of neutron stars, at least with ref- 
erence to one of its possible consequences, namely Ekman pumping. This 
is a mechanism that can considerably shorten the timescale needed for the 
redistribution of angular momentum (in comparison with viscous diffusion 
characterized by the timescale given by r visc ps Rljv*. where v* is the typical 
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kinetic viscosity coefficient and i£* is the relevant stellar radial length scale) 
and thus the damping of differential rotation in cases for which (as will be the 
case in a typical pulsar) the star is rotating fast enough for the correspond- 
ing rotation timescale r rot = 2n /VL to be short compared with r visc . In such 
circumstances, "Ekman pumping" will supplement the very slow diffusive 
transport by more rapid convective transport propelled by centrifugal buoy- 
ancy forces. The ensuing "Ekman timescale" r E for the effective damping 
of differential rotation in such cases will be given roughly by the geometric 
mean of the pure diffusion and rotation timescales, i.e. r E « yY rot T v i sc . 

While it has been recognized that either Ekman pumping or magnetic 
coupling is in general efficient to bring into corotation the core plasma with 
the crust |TJ, it is expected that Ekman pumping is quite inefficient (see 
e.g. 0) for the uncharged crust neutron superfluid that is believed (see 
e.g. P) to permeate the lower layers of the crust in the density range from 
10 11 to about 10 14 gm/cm 3 . This means that the convectively accelerated 
Ekman timescale, r E R^^2ttJuSI, is too long to prevent the development 
of significant differential rotation. The negligibility, in such cases, of Ekman 
pumping is attributable to the effective negligibility of viscosity, but should 
not be construed as implying the negligibility of centrifugal buoyancy forces. 
In previous discussions of such scenarios - and in particular of the simplified 
strictly stationary limit in which the effective viscosity is neglected, so that 
no possibility of Ekman pumping can arise at all - the role of centrifugal 
buoyancy forces has been rather generally overlooked. The upshot of the 
present investigation of stationary differentially rotating configurations is to 
show that in such cases the general neglect of the centrifugal buoyancy effect 
is quite unjustified, and that on the contrary this effect is potentially capable 
by itself of providing the dominant contribution to the crust stresses that are 
ultimately released in "glitches" . 

2 Glitches driven by the spheroidality mech- 
anism. 

In the first years after the problem of accounting for neutron star glitches was 
posed, attention was concentrated on what is describable as the "spheroidal- 
ity mechanism" ||. This mechanism depends on the supposition that the 
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solidity forces will not be strong enough to allow the stellar equilibrium con- 
figuration to differ very much from a perfectly fluid equilibrium state, which 
would be spherical in the absence of rotation, but which will actually have the 
form of an oblate spheroid with ellipticity proportional to Q 2 . The moment 
of inertia, defined as the ratio I = J/Q, will be given for a slowly rotating 
fluid by an expression of the form 

/ = / (l + fi 2 /U 2 ), (1) 

where I Q is its spherical limit value and is a constant characterising the 
rather high angular frequency needed for relative deviations from spherical 
symmetry to be of order unity. A more accurate formula involving higher 
order corrections would be needed for a star with angular velocity near the 
critical value Q 2 ~ Q 2 , but the cases in which glitches have been observed so 
far are all characterised by 

fi 2 « V 2 . (2) 

For a perfectly fluid star model, a continuous angular momentum variation 
AJ < would bring about a corresponding momentum of inertia variation 
AI that would be given by 

AI Q 2 Att 

- S32 njir <0 - < 3 » 

Due to the solidity of the crust, which tends to preserve the more highly 
elliptic initial configuration, the actual change in the moment of inertia will 
fall short of what is predicted by this formula, but at some stage the strain 
will build up to the point at which the solid structure will break down (see 
Fig. 1) . It is predicted that there will then be a "crustquake" in which the 
solid structure suddenly changes towards what the perfect fluid structure 
would have been, thereby changing the moment of inertia by an amount 

51 = eAI, (4) 

where e is an efficiency factor that should presumably lie somewhere in the 
range 

0<£<1. (5) 

Since the amount of angular momentum loss during the very short duration of 
the glitch will be negligible, the corresponding discontinuous angular velocity 
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change will be given by 

IT = ~T' (6) 

Its value will therefore be expressible in terms of the order of unity efficiency 
factor e by 

5n = -2^^Afi , (7) 

in which it is to be recalled that AQ denotes the continuous (negative) change 
in angular velocity since the preceding glitch. This mechanism must presum- 
ably operate, and may account for some observed glitches, but it soon became 
clear [|[ that even if this mechanism is maximally efficient, with 

e ~ 1 , (8) 

the magnitude predicted by (0) is much too low for such a mechanism to 
be able to account for the comparatively large glitches that are frequently 
observed in cases such as that of the Vela pulsar. 



3 Glitches driven by differential rotation. 

Soon after the empirical discovery of glitches too large to be accounted for 
by the "spheroidality" mechanism, it came to be recognised by theorists p] 
that a plausible explanation involved the superfluid property of the deeper 
layers of sufficiently cool neutron stars. This property makes it possible to 
conceive that an interior neutron superfluid layer with moment of inertia, J n 
say, can rotate with an angular velocity, Q n say, that may differ from the 
externally observable angular velocity Q that characterises the part of the 
star that corotates with the crust, with its own moment of inertia 

I C = I-I Q . (9) 

In such it can be supposed that when an external braking mechanism 

causes the corotating crust component to undergo an angular velocity change 
AQ, the angular velocity Q n of the independently rotating neutron superfluid 
layer may in the short run be unaffected, with negligible variation expressible 
by 

AQ n = 0, (10) 
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but that, when the ensuing angular velocity difference between the corotating 
crust component and the neutron superffuid layer exceeds some critical value 
there will be a discontinuous adjusment whereby this angular velocity differ- 
ence is reduced by some process involving a transfer of angular momentum 
between the two components. Such a process will evidently entail a negative 
adjustment 5Q n of the angular velocity of the neutron superffuid layer and an 
accompanying positive adjustment 8Q of the (observable) angular velocity of 
the corotating crust component, whereby the latter increases its angular mo- 
mentum by an amount I C SQ that is equal to the amount —I n SQ n that is lost 
by the neutron superffuid component, so that the total angular momentum 
change during the discontinuous 'glitch' process is zero, i.e. 

i c sn + ijn n = o . (ii) 

If this adjustment process were a hundred per cent efficient, the net variation 
AQ + Sfl of the corotating crust angular velocity would be exactly matched 
by the net neutron superffuid angular velocity variation, which by (|T0| ) will 
be simply given by 5Q n , so that one would have 

5tt n = e(Att + 50.) , (12) 

with e ~ 1. In practice one would expect that there would typically be an 
incomplete adjustment, still expressible by a relation of the form (VI), but 



with an efficiency factor e having some lower value in the range (|5|). By 
substituting (|IID in ( |12|) it can be seen that the observable glitch magnitude 
will be given by 

m = T^r;r (13) 

and hence, by (|9|), that for an efficiency factor e with any value in the range 
(H) the glitch magnitude will satisfy the inequality 

5n>-£jAn. (14) 

By comparing ([H]) with (|7|), it can be seen that, for a given assumed 
value of the efficiency factor e, the differential rotation adjustment mecha- 
nism characterised by (|T4]) can give rise to a much larger glitch magnitude 
5Q than is possible by the spheroidality adjusment mechanism characterised 
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by (|7|), because the factor I n / 1 in (TfiJ) can be of order unity, whereas the 
corresponding factor in (0), namely 2Q 2 /Q% is very small compared to unity 
in even the most rapidly rotating pulsars. Thus, unlike the spheroidality 
mechanism, mechanisms involving angular momentum transfer between dif- 
ferentially rotating components can plausibly be considered as candidates for 
explaining the frequent large glitches observed in the Vela pulsar. 



4 Glitch mechanisms due to the vortices. 

In the context of a glitch due to differential rotation, the question that arises 
is what physical mechanism can increase the effective coupling between the 
superfluid component and the crust, in order to generate a transfer of angular 
momentum. 

The explanations that exist in the literature are based on an important 
property of a superfluid neutron star, which we have not yet mentioned in this 
article: the existence of an array of vortex lines in the rotating neutron super- 
fluid component, each vortex carrying a quantum of vorticity k = h/(2m n ) 
(where m n is the neutron mass). The vortex number density (per unit area) 
n v is directly related to the superfluid angular velocity Q n by the expression 

n v = (15) 

K 

(for uniform rotation). 

The kind of angular momentum transfer mechanism that has for many 
years been generally considered to offer the most likely explanation for large 
glitches is based on the supposition that these vortices will be "pinned" in the 
sense of being effectively anchored in the lower crust, either by pinning in the 
strict sense or by a sufficiently strong friction force 0. The braking of the 
crust will thus have the effect of slowing down the vortices relatively to the 
underlying superfluid, thereby giving rise to a Magnus force tending to move 
them out through the superfluid layer and thus slow it down as well. However 
this tendency to move out will be thwarted by the same anchoring effect that 
gave rise to it in the first place. This conflict will cause the pinning forces 
to build up to a critical point at which there will be a breakdown bringing 
about a discontinuous readjustment of the kind described by the analysis of 
the preceding section, and in particular by the formula (|l4|). 
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The breakdown can occur in two different manners: 

(a) There can be a sudden unpinning of many vortices, due to the break- 
ing of the pinning bonds 0,0. 

(b ) Another possibility is that the crust lattice breaks before vortex lines 
can unpin from it, as suggested in || and studied in detail by Ruderman [|]. 

Finally, we would like to mention another interesting glitch mechanism 
due to Link and Epstein which may be relevant for the present work: 

(c) their thermally driven glitch mechanism is based on the so-called 
vortex creep model [0], in which the coupling between the vortices and the 
crust is strongly temperature dependent. A sudden local increase of the inner 
crust temperature, such as may be due to a crustquake, can then be shown 
to induce a glitch. 

It must be emphasized that all these three mechanisms, even if corre- 
sponding to some breaking of the crust as in the scenarios (b) and (c), are 
very different from the mechanism of section 2, in the sense that they all are 
in the context of a two-component star, with the neutron superfluid rotating 
faster than the crust and thus acting as a reservoir of angular momentum. 
In the following sections, we will consider a mechanism which is not based 
on the presence of vortices, but still in the context of differential rotation. 

Finally, let us mention the question of how big is I n compared with /, 
in other words how much of the neutron fluid is effectively free to rotate 
independently of the rest? In the unpinned part the vortices can move out 
freely so as to establish corotation, so I a may be relatively small [111, repre- 
senting the moment of inertia just of the small fraction of the neutron fluid 
that interpenetrates the deeper layers of the solid crust where pinning is ex- 
pected to be most effective. However effective pinning may not be confined 
to the solid crust: it may also be achieved by forces exerted by quantised 
magnetic field lines (resulting from superfluidity of the protons) in the layers 
below the crust, in which case the relevant value of I n might be much larger 
[ |T2] ] . Another question (which applies also to the less important spheroidality 
mechanism discussed above) is that of the absolute values of the discontin- 
uous changes. The foregoing reasonning is concerned just with the ratio of 
5Q to — AQ but does not tackle the harder problem of their absolute values. 
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5 Potential importance of the centrifugal buoy- 
ancy mechanism. 

So far we have only been summarising what has long well known to workers 
in this field. We now come to what seems to us to be an important point that 
has been overlooked, which is that independently of vortex pinning there is 
another, comparably powerful mechanism, that can also cause discontinuous 
angular momentum transfer to a solid crust from an independently rotating 
superfluid layer. This mechanism does not depend on superfluidity in the 
strict sense but merely requires perfect fluidity in the sense that the effective 
viscosity should be low enough for the slowdown of the neutron fluid to lag 
behind the slowdown (due to its coupling with the radiating magnetosphere) 
of the solid outer layers. The point is that if the outer layers were also effec- 
tively fluid, there would be a convective readjustment, in which annular rings 
of fluid would change their relative positions, each retaining its separate an- 
gular momentum, in such a way that those with less angular momentum per 
unit mass, and thus with less "centrifugal buoyancy" would move towards 
the axis while those with more would move out so as to establish a state 
of equilibrium in which, provided the pressure depends only on the density, 
the angular velocity would decrease outwards as a function just of cylindrical 
radius, in accordance with the well known Taylor- Proudman theorem (see, 



e.g., [p!3[ ). The effect of crust solidity will be to temporarily postpone such 
readjustments, by the development of the anisotropic stresses needed to bal- 
ance the centrifugal buoyancy forces. However when such stresses have built 
up to the critical point at which the solid structure breaks down, the pent 
up centrifugal buoyancy forces will produce a "starquake" in which the con- 
vective readjustment that would have ocurred continuously in the fluid case, 
is finally achieved in a discontinuous transition. 

It is to be noticed that in contrast with the vortex pinning effect (in the 
following, for easier comparison, we will have in mind the scenario (b) of Sec- 
tion 4), which tends to pull the more slowly rotating crust material outwards 
from the axis towards the equator (see Fig. 2), the effect of the centrifugal 
buoyancy deficit in the crust is to pull the crust material inwards towards 
the axis of the star, where it will finally be subducted into the fluid interior 
(see Fig. 3). Although the centrifugal buoyancy effect produces convective 
circulation in just the opposite direction to that produced by vortex pinning 
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(which if it were strong enough would lead to subduction at the equator 
rather than the axis ||14|| ) its effect on the angular momentum distribution 
would be similar, i.e. the net effect of a centrifugal buoyancy crustquake will 
be a discontinuous transfer angular momentum to the crust from the more 
rapidly rotating fluid layer. This means that the crude quantitative estimate 
given by equation fll4|) is applicable just as well to the effect of a centrifugal 
buoyancy crustquake as to a vortex pinning crustquake. 

The main point we want to emphasise is that whereas vortex pinning 
may indeed be the main driving force for the build up of the stress that is 
relaxed in crustquakes, the extent to which it really is depends on detailed 
considerations about the strength of vortex pinning. On the other hand 
the opposing centrifugal buoyancy mechanism will always function whenever 
there is differential rotation. It will be seen in the next section that when 
it is fully effective the oppositely directed pinning mechanism will be strong 
enough to overwhelm (i.e. to more than cancel) the buoyancy mechanism, 
but the latter mechanism is more robust in the sense that it will always make 
a significant contribution. 

Our tentative conclusion - which we are proposing as a subject for debate 
and further investigation - is that the hitherto neglected centrifugal buoyancy 
effect may be the dominant cause of the crustquakes that are observed as 
pulsar glitches, while vortex pinning crustquakes, if they occur at all, are 
relatively rare. This does not mean that vortex pinning is unimportant for the 
phenomenon, because it is likely to be what determines the magnitude of the 
relevant moment of inertia contribution J n in the estimate (0) for the ratio 
of 5Q to —AQ. However what it means is that the vortex pinning stresses are 
not what is immediately responsible for the discontinuous breakdown, and 
hence not what is of dominant relevance for estimating the absolute values 
of AQ at which it is likely to occur. 

6 The working of the centrifugal buoyancy 
deficit mechanism. 

An accurate treatment of neutron star would of course require a general 
relativistic analysis ]15, 16], but as a first step towards the estimation of the 



stress forces needed to maintain equilibrium where the crust constituent is 
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interpenetrated by an independently rotating fluid constituent, it will suffice 
for our present purpose to work in a Newtonian framework, using a highly 
idealised two-constituent model in which the corotating crust component 
(including the protons and electrons, as well as a fraction of the neutrons that 
is bound into atomic type nuclei) and the neutron superfluid are considered 
as independent material media having respective mass densities 

p c = mn c , p n = mn n (16) 

and spatial velocity components v£ and v£ (i = 1, 2, 3) where m is the proton 
mass and n c and n n are the corresponding baryon number densities. For an 



approximate description of the kind of scenario envisaged by Alpar et al [17 
in which the rigidly corotating constituent consists not just of the crust lattice 
but also of the proton superfluid in the core which will be locked to the crust 
by electromagnetic interactions we adopt a simplified treatment in which it is 
postulated that the dynamics is governed by Euler type equations of motion 
of the familiar form 

Pc {dy c + vi^jv*) = -v l p c - Pc v> + (17) 

Pn (d < + vlVjV*) = -V*P n - Pn VV + fi (18) 

using d to denote partial differentiation with respect to Newtonian time, 
where <fi is the Newtonian gravitational potential, and where P c , P n and 
/ c l , /jj respectively denote the relevant pressure scalars and force density 
vectors. In a lowest order approximation in which both components can be 
considered to obey barotropic equations of state giving their energy densities 
e c and e n as functions respectively of n c and of n n , they will be characterised 
by corresponding chemical potentials 

d,E c d,£ n /tr>\ 

^ c= d^> ^ n= dn~' (l9) 
from which the associated pressure contributions can be evaluated as 

P c = Pc^c -e c , P n = p n n n - e Q . (20) 

This implies that the required gradient terms will be given by 

V l P c = n c V> c , V*P n = n n V> n . (21) 
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(It is to be remarked that in a more detailed analysis the baryon chemical po- 
tential /j, c in the component corotating with the crust would be interpretable 
as the sum of proton and electron contributions, \i c = /i p + /i e .) 

Although adequate for the fluid constituent, a purely barotropic descrip- 
tion will not be sufficiently accurate for the crust constituent in which we 
want to allow for the effects of solidity. The usual way to do this is to replace 
the isotropic pressure gradient term VjP c by a stress gradient term of the 
form VjT,?i where is the total stress tensor. It will be convenient for our 
purpose to decompose the latter in the form 

T C { = P C S{- si (22) 

where the extra anisotropic stress contribution s\ is a correction term that 
will be small compared with the dominant isotropic contribution P c 5j. This 
means that while / n l is to be interpreted as the interaction force density, if 
any, exerted on the neutron superfluid component by effects such as vortex 
pinning, on the other hand the term / c * in ( [T7D will consist, not just of 
the equal and opposite interaction term — / n J but also of an extra correction 
term f s l due to the anisotropic stress correction representing the effect of the 
solidity property, i.e, we shall have 

fi = fi-fi, /.' = -V,. (23) 

The anisotropic stress contribution s\ and the associated force density / S J 
might also include an allowance for magnetic effects, such as are ultimately 
responsible for the external braking mechanism and for locking the proton 
superfluid in the core to the outer crust lattice. However for the equilibrium 
of the strictly stationary states with which we shall be concerned here such 
magnetic effects are not important, so it may be considered that the stress 
force density f s l arises just from the Coulomb lattice rigidity in the crust, 
and that it vanishes in the high density core. 

Let us now restrict our attention to configurations that are stationary, 
so that the terms acted on by d will vanish, and let us suppose the motion 
consists just of a circular motion about the x 3 axis, so that each comoving 
particle moves with a fixed value of the cylindrical radius w = (x l2 + x 2 " 2 ) 1 ^ 2 . 
This means that the velocity gradient terms in the equations of motion will 
be given by 

viVjvl = -\tf c VW , vlVX = -\^ n VW , (24) 
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where f2 c is the local angular velocity of the crust constituent and Q n is the 
local angular velocity of the superfluid constituent. Under these conditions 
the Euler equations (|T7p and QI%D can be rewritten in the form 



'^W - V*(</> + m"Vc) = p-'if: - //) , (25) 



2 

and 

2 



X fi n 2 VV 2 - V*(0 + m-Vn) = -ftTVn* • (26) 



If vortex pinning were effective, it would contribute to / n l the force density 
needed to counteract Joukowsky-Magnus type lift force density fj that would 
be exerted on the vortices by the Magnus effect, which would be given by 

fj = Pn (n n -n c )n n v l w 2 , (27) 

but in the absence of vortex pinning or other coupling forces, the right hand 
side of (|26|) will simply vanish, in which case it can be seen that the fluid will 
satisfy the Taylor-Proudman condition, meaning that its angular velocity 
Q n and also the combination mcj) + p n must vary as a function only of the 
cylindrical radius w. 

Since the interaction force density / n l will cancel out of the linear com- 
bination of fl25D and (E?6| ) obtained from the direct sum of ( |17D and (|18D , it 
follows that this combination will take a simple form that is conveniently 
expressible - independently of whether vortex pinning is actually effective or 
not - in terms of the "would-be" Joukowsky force density ( p7|) as 



v*p + P (w - |ft 2 w 2 ) = f* + /; - i Pn (fi n - q c ) 2 w 2 , (28) 

in which the total pressure P and mass density pressure p are defined in the 
obvious way as 

P = P c + P D , p = p c + p n . (29) 

In a systematic calculation by successive approximations, the first stage 
would be to obtain a zeroth order solution of the stellar equilibrium problem 
in which the (first order) crust rigidity and differential rotation contributions 



on the right hand side of (|28|) would simply be neglected. What we are 
interested in here is the next stage, which involves the first order equation 
(from which the zeroth order part has cancelled out) that is obtainable by 
taking the difference of (p5|) and (E? 
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Before going ahead it is necessary to stress that, since only weak inter- 
actions are involved, it cannot be taken for granted that the relevant nu- 
clear transitions involved in the "neutron drip" process whereby matter is 
transferred between the ionic crust material and the interpenetrating neu- 
tron superfluid will be very rapid compared with the "secular evolution" 
timescales on which the state under consideration is significantly modified. 
If the "neutron drip" process were sufficiently rapid one would obtain not 
just mechanical equilibrium, such as expressed by equations ( p5|) and (|26|), 
but also thermodynamical equilibrium in the rest frame of the crust, in the 
sense that the energy per baryon of the "normal" matter corotating with 
the crust, which is just /i c , would be the same as the energy per baryon of 
the neutron fluid with respect to the crust corotating frame, which has the 
value /i n + \m{Vt n — Q c ) 2 zu 2 . In practice however, due to the slowness of the 



relevant nuclear transitions [TE , it is necessary to allow for the possibility of 
a finite deviation, 

Afj, = fj, c - fj, n - -m(fi n - Q c ) 2 w 2 , (30) 

from exact thermodynamic equilibrium. Estimates of the likely values for 
such a chemical potential excess due to the simple spheroidality adjustment 
mechanism, discussed above in Section 2, have been provided by the recent 
work of Reisenegger [PJ. Significantly larger values are likely to arise from 



the differential rotation mechanisms considered here due to the resulting 
tendency for the crust constituent to be convected relative to the neutron 
fluid constituent. 

Including allowance for the possibility of a neutron drip delay contribution 

fi = n c V*(A/i) , (31) 

representing the force density due to the chemical potential excess (ffflp if 
any, the solid stress force density / s l ultimately responsible for the glitches in 
which we are interested can be seen to be given by the first order equation 
obtained by subtracting (EH) from (Bop, which will be expressible in the form 



fi = fi + f fi + fi • (32) 

Pn 

The final term in the above equation is what can be interpreted as the extra 
force needed to compensate for the buoyancy deficit of the crust due to its 
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lack of rotation velocity relative to the neutron superfluid, and is given by 

/„* = Pc (fi n " fi ) (V* ( U7 2 n n ) - W 2 V ; fi c ) . (33) 



7 Estimation of the centrifugal buoyancy deficit 
force density. 

The solidity property of the crust implies that, in a stationary state, its 
rotation must be rigid, i.e. 

tt c = n , V'ft = , (34) 

where Q is a uniform angular velocity value (the one that is actually observ- 
able from outside), so the formula ([33]) for the buoyancy deficit force density 
can be immediately simplified to the form 

/ b l = Pc (fi n - n)v l (w 2 n n ) . (35) 

If the superfluid were macroscopically irrotational, i.e. if there were no vor- 
tices present, then w 2 Q n would have a uniform value so the right hand side 
of (|35D would also vanish, i.e. the effective buoyancy deficit force density / b J 
would be zero. 

What we actually anticipate in the context of the pulsar slowdown prob- 
lem is that Q n will be approximately uniform (representing rigid rather than 
irrotational motion) with a value equal to that of the crust component at 
a rather earlier stage, perhaps just after the previous glitch, and that the 
velocity difference will therefore be small compared with the total angular 
velocity 

|fi n -fi|<|0|. (36) 

Thus, by neglecting corrections of quadratic order in this velocity difference, 
we see that Q3"5]) can be conveniently approximated by the simpler formula 

ti~ Pc (n n -n)n n v\w 2 ), (37) 

which will be accurate to linear order in the difference Q n — Q. It is to be 
remarked that, to the same order of accuracy, the neutron drip delay force 
contribution ([H]) will be given by the approximation 

f:~n c V\p c -p n ). (38) 
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It is to be observed that the formula ( |37[ ) for the buoyancy deficit force 
density closely ressembles the Joukowsky formula (p7|) for the lift force den- 
sity /j that would be exerted on the vortices by the Magnus effect if they 
are pinned to the crust: this Joukowsky force density is evidently related to 
the buoyancy deficit force density by the simple proportionality relation 

fj * ^/b • (39) 

Pc 

It follows that, in terms of the effective (centrifugally adjusted) gravitational 
potential 

= - , (40) 

the basic stellar equilibrium equation (|28|) will reduce to the form 

v*p + ~ n + /; , (4i) 

in which the zeroth order terms are grouped on the left and the first order 
terms are on the right (while the final second order term on the right of 
( P5[) has been neglected). Since the left hand side consists just of the small 
difference left over after the approximate cancellation of the dominant zeroth 
order terms, this equation does not provide any utilisable information about 
the solid force density / s l in which we are interested: on the contrary, after 
/ s l has been evaluated by other means, (|4lD can be used to calculate the 
corresponding first order adjustments to the zeroth order pressure and density 
distributions. 

The equation that does supply the relevant information about the solid 
stress force density / s * in which we are interested is the first order equilibrium 
condition (|32|), whose terms can be instructively regrouped in the form 

/; - fi = /b + f fl , (42) 

Pn 

in which it can be seen from (|37D that the right hand side will always be 
approximately proportional to the first order difference Q a — Q whether or 
not pinning is effective. (This shows incidentally that differential rotation 
would be impossible if both the rigidity force / s l and the chemical delay 
contribution / x * were negligible.) 



16 



In particular, the relation (pE2| ) shows, by ( |39| ) that in the pinned case, i.e. 
when the superfluid is submitted to a force density 



f: * -fi ■ (43) 

on the crust, the stress force density / s l necessary for equilibrium will be 
given by the simple formula 

/■* = -n + fi , (44) 

in which the first term on the right is just the Joukowsky- Magnus contribu- 
tion, as is assumed in the conventional presentation of the vortex pinning 
theory of pulsar glitches. 

The formula (fP|) is potentially misleading in that it gives the false im- 
pression that if the pinning were ineffective, so that instead of being given 
by fl43D the force exerted on the superfluid by the crust were simply zero, 

fn = , (45) 

then the first term on the right of the stress force density formula would 
similarly disappear, whereas in fact substitution of ( [45| ) in (|32| ) leads to the 
replacement of (f0]) by the formula 

/; = k + fi , (46) 

in which, instead of the Joukowsky- Magnus contribution — /j , the right hand 
side is now given by the oppositely directed buoyancy deficit force contribu- 
tion / b J . 

Our reasonning so far does not make it obvious whether or not the crust 
will develop a sufficiently non-uniform chemical potential excess A/i to pro- 
vide a significant chemical excess force / x \ If it is a good approximation to 
suppose that chemical excess force in the crust vanishes, 

fl = , (47) 

(as seems to have been implicitly asumed in most previous works but which 
needs to be confirmed or infirmed quantatively) then, in the case where 



pinning would not be effective (as has been advocated by Jones |IT] contrarily 
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to earlier works), it follows from ( pE6| ) that there will still be a solid stress 
force density given by 

* /b , (48) 
in which the centrifugal buoyancy deficit force density on the right is given 
by equation (PT|). This formula can be seen to differ from the (alternative) 
well known formula - for the stress due to pinning -, deduced from fl44|) by 



the same assumption (|47|), 

ft * ~fj , (49) 
with the Joukowsky- Magnus term on the right hand side given by (P7|), only 
by having the opposite sign and by having a proportionality factor given by 
the density p c of the corotating crust component instead of the density p n of 
the differentially rotating neutron superfluid component. 



8 Discussion and conclusions. 

In the lower crust region that seems most likely to be relevant for the ex- 
planation of the large glitches observed in the Vela pulsar one would expect 
the corotating constituent to be characterised by a density p c (attributable 
mainly to protons and bound neutrons in the atomic type ions forming a solid 
lattice) having a range of values that is roughly comparable with that of the 
corresponding neutron superfluid density p n (quantitatively round about 10 13 
g/cm 3 ). Thus although they are of opposite sign (tending to push the crust 
material outward in the case ( |43|) of vortex pinning, but to push it inwards in 
the case ( [45]) for which pinning is absent) the alternative formulae ( fj8|) and 
( fl9|) both predict the same rough order of magnitude for the stress induced 
on the crust by the existence of a difference between the angular velocity Q n 
of the neutron superfluid constituent and the (externally observable) angular 
velocity Q characterising the crust. 

The implication is that, as a candidate for explaining the large magnitude 
of the discontinuous changes 5Q that are commonly observed in a pulsar such 
as Vela, the previously overlooked buoyancy deficit mechanism characterised 
by the formula fl48|) , i.e. 

f:^p c (Q Q -Q)Q n V^ 2 , (50) 

(pushing outward along the cylindrical radial direction) seems at first sight 
to be just as promising as the more thoroughly investigated vortex pinning 
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mechanism, which, if the chemical contribution / x * were unimportant, would 
be given according to (|49l) by 

f:~-p n (n n -n)n n v l m 2 , (51) 

(pushing inward along the cylindrical radial direction). In order to obtain 
definitive conclusions it is clear however that much more work on both kinds 
of mechanism will be needed. In particular it will be necessary to pay more 
attention than hitherto to the role of the chemical excess force fl3"T|). 

The present situation can be summarised by the statement that the large 
magnitude of the observed glitches in Vela provides strong evidence for the 
existence of angular velocity differences - and hence for the existence of 
superfluidity - in the pulsar interior, but that it is premature to claim it also 
provides strong evidence for vortex pinning because stresses of comparable 
magnitude could be produced in the absence of pinning by the centrifugal 
buoyancy deficit mechanism. 
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Figure 1: Qualitative sketch indicating direction of force expected to act on 
(magnetically slowed down) crust due to spheroidality mechanism, in absence 
of differential rotation. Solid lines indicate outer and inner boundaries of 
crust. Vertical shading indicates the alignement of the vortices in the region 
occupied by neutron superfluid, which is not confined to the core but inter- 
penetrates the greater part of the solid cust as well. Note that the vortices 
represented here are not physically relevant in this particular mechanism. 




Figure 2: Qualitative sketch indicating direction of force expected to act on 
(magnetically slowed down) down on crust due to vortex pinning mechanism, 
if it is effective, when the (interpenetrating) neutron superfluid retains a 
higher rotation rate. 
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Figure 3: Qualitative sketch indicating direction of force expected to act on 
(magnetically slowed down) crust, even if vortex pinning is ineffective, due 
to the centrifugal buoyancy mechanism when the (interpenetrating) neutron 
superfluid retains a higher rotation rate. 
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